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ABSTRACT 



ISO has provided us with a new perspective on gas 
temperatures in photodissociation regions through 
measurements of hue emission from rotationahy- 
excited levels of H2. The H2 rotational level pop- 
ulations provide a thermal probe, showing that gas 



temperatures Tg, 



500 — 1000 K prevail in a por- 



tion of the PDR where significant H2 is present. Such 
high gas temperatures were unexpected. Theoretical 
models for the S140 PDR are presented. Possible 
mechanisms for heating the gas to such high temper- 
atures are discussed. 



Key words: ISO; infrared astronomy; molecular hy- 
drogen. 



1. INTRODUCTION 



Photodissociation fronts play an important part in 
the global energetics of star-forming galaxies, as 
an appreciable fraction of the energy radiated by 
newly-formed stars is reprocessed in photodissocia- 
tion fronts, resulting in H2 and fine structure line 
emission. If we want to understand the emission spec- 
tra of star-forming galaxies, good theoretical models 
of photodissociation fronts are required. 

ISO observations of line emission from rotationally- 
excited levels of H2 have provided unequivocal evi- 
dence for gas temperatures in the 500 — 1000 K range 
in a portion of the photodissociation region (PDR) 
where the H2 fraction is appreciable. These temper- 
atures were higher than expected from current mod- 
els of the heating and cooling processes in PDRs, and 
therefore require reconsideration of the physics of the 
gas and dust in PDRs. 

In §p] we review the basic structure of photodissocia- 
tion fronts, and in §^ we discuss how H2 can be used 
as a thermometer to indicate the gas temperatures 
in PDRs. In ffl we use the S140 PDR as an exam- 
ple. Theoretical models for the heating and cooling 



in PDRs are discussed in 
some of the uncertainties in 
We summarize in §7] 



, and in ^ we stress 
le modelling of PDRs. 



2. WHAT IS A PHOTODISSOCIATION 
FRONT? 



The term "photodissociation front" refers to the in- 
terface separating a region which is predominantly 
molecular, and a region where the ultraviolet en- 
ergy density is sufficiently high so that the molec- 
ular fraction is <c; 1. Young O and B stars located 
near their natal (only partially-disrupted) molecular 
clouds produce conspicuous (i.e., high surface bright- 
ness) photodissociation fronts. A number of reviews 
of PDRs have appeared recently ( Hollenbach fc Tie - 



lens 1997, Sternberg et al 199g |Walmslev jgTTR 



mpreJiensive article by hLollenbacJi fc 



rielens (1998)| . 

A photodissociation front is not a mathematical sur- 
face, but rather an extended "photodissociation re- 
gion" (PDR). Figure H is a cartoon illustrating the 
different zones within the PDR. 

The structure of the PDR is determined primarily by 
the attenuation of the far-ultraviolet (6 - 13.6 eV) ra- 
diation field, as one moves from the ionization front 
into the PDR. The H2 abundance results from a bal- 
ance between formation of H2 on dust grains and pho- 
todissociation of H2 by 912 < A < lllOA photons. 
The photodiss ociation rate is determin ed by both H2 
self-shielding (praine fc Bertoldi 1996 ) as well as at- 
tenuation by dust. Because ot this self-shielding, the 
H/H2 transition occurs closer to the ionization front 
than the C+/C/CO or O/O2 transitions. Atoms with 
ionization potentials lower than H (e.g., C) are pho- 
toionized in the PDR, whereas species with ionization 
potentials greater than H (e.g. O) are neutral in the 
PDR. 

Ionization fronts and dissociation fronts are, in gen- 
eral, time-dependent structures ( Bertoldi fc Draine 
1996). Under many circumstances, however, the 
structure can be approximated as being station- 
ary in a frame moving with the dissociation front. 
Furthermore, if the dissociation front is advancing 
into the molecular gas with a propagation speed 



^DF ^ 0.3 kms (Bertoldi fc Draine 1996), advec- 



tion terms can be neglected in the equations govern- 
ing the thermal and chemical conditions: the tem- 
perature is effectively determined by a local balance 
between heating and cooling, and abundances by a lo- 
cal balance between formation and destruction. For 



most PDRs this inequality is expected to be sat- 
isfi ed (pertoldi fc Draine 1996 , [Storzer fc Hollen^ 
ba( :h l&98| j, altiiougii time dependen ce may be im- 



portant tor some planetary nebulae ( Hollenbach 
N 4ta 1995| ). 



HII 



Photodissociation Region 




Figure 1. The structure of a PDR produced by radi- 
ation from a star which also emits enough ionizing 
photons to produce an HII region. The ionization 
front separates the HII region from the PDR. Typi- 
cal values of the effective far- ultraviolet optical depth, 
and H nucleon column density are shown. 



3. H2 AS A THERMAL PROBE 

The (v, J) royibrational excited states of H2 are pop- 
ulated by inelastic collisions 

H2(^i, Ji) + X ^ H2(^2, J2)^X±KE ; 

by UV pumping: 

H2(0,Ji) + /i^ ^ H*(^i,Ji±l) 

^ H2(^2,Ji±0,2) + /i^(UV) 
^ H2(^3,Ji ±0,2,4) + /ii/(IR) 

^ H2(0,Ji±0,2,4,...) ; 

and by formation on grains: 

grain + H + H ^ grain + H2 (v, J) + KE . 

At the densities nu ^ 10^ cm~^ of bright PDRs, col- 
lisions maintain the (v=0, J) leyels of H2 in approxi- 
mate thermal equilibrium for J ^ 5. Therefore mea- 
surements with ISO of the quadrupole emission line 
intensities from leyels J ^ 5 proyide a good indicator 
of the gas temperature, while the emission from the 
higher J leyels is sensitiye to both the gas temper- 
ature and density. We do not expect the leyel pop- 
ulations to be characterized by a single "excitation 



temperature" , but the "best-fit" excitation tempera- 
ture characterizing a range of J yalues (e.g., J=3-5) 
should indicate the approximate gas temperature in 
the part of the PDR where these leyels are predomi- 
nantly excited. 

The populations of yery high J leyels (e.g., J = 
15) are potentially affected by UV pumping (the 
quadrupole decay cascade following a UV pump in- 
jects a small amount of H2 into high J states) and 
formation on dust grains (some fraction of the newly- 
formed H2 may be in high J states). Therefore ex- 
citation temperatures characterizing J ^ 10 may not 
reflect gas temperatures. 



4. THE S140 PDR 



ISO obtained spectra of the photodissociation region 
where the S140 HII region abuts the L12 Q2/L1204 
molecular cloud ( Timmermann et al. 199q ); see Fig- 
ure 0. 
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Figure 2. The SI 40 PDR, where radiation from 
the B0.5V star HD 211880 impinges on a dense 
clump in the L1202/L1204 molecular cloud. The 
ionization front is shown as a shaded region. The 
grey scale shows [CI]610/j.in emission (White & Pad- 
man 1991). Contours show CO (3-2) emission 
(Hayashi et al. 1987). The SWS apertures are shown 
as rectangles positioned just interior to the ionization 
front, where the photodissociation front is expected to 
be located, viewed approximately edge- on. The excit- 
ing star is located ~7 arcmin SW of the aperture. 
Figure taken from Timmermann et al.(1996). 

I I 

Timmermann et al. (1996) reported measurements 
of line emission out of (0,J) leyels with J as large 
as 9 (see Figure ^). The inferred column densities 
for 3 < J < 7 are characterized by an excitation 
temperature Texc ^ 500 K; the emission from J = 8 
and 9 suggest a higher temperature. How much of 
the PDR is characterized bv such high temperatures? 
To address this question, Timmermann et al. (1996)| 



considered an ad-hoc temperature profile, shown in 
Figure ^. 
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[Sill] 35 /im and [Fell] 26 fivn agree with observations if 
Si/H=1.5 X 10-^ and Fe/H=5 x 10"^ The principal 
discrepancy is [OI]63/im, which the model predicts 
to be ^ 5 times stronger than observed. 
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Figure 4- Ad-hoc temperature profile adopted for SI 40 
PDR model. With this assumed temperature pro- 
file, the H2 abundances are computed by requiring a 
balance between H2 formation and photodissociation. 
Taken from Timmermann et al.(1996). 



Figure 3. Observed (filled symbols) and model (open 
symbols) column densities in rovibrational levels of 
H2 along the line- of- sight. Observed intensities have 
been corrected for foreground extinction by dust with 
0.3 mag of extinction at K. The PDR is assumed 
to be inclined with cosO = 0.1. The broken line 
is a T = 500 K thermal distribution, indicating 
that a significant amount of H2 is located in regions 
where Tgas ^ 500 K. Taken from Timmermann et al. 
(1996). 

With this ad-hoc temperature profile, the chemi- 
cal abundances are computed by requiring steady- 
state balance between formation and destruction, us- 
ing standard assumptions for the formation rate of 
H2 on grains, and a detailed treatment of H2 self- 
shielding. H2 level populations were computed for 
a plane-parallel PDR, including both UV pumping 



an d collisiona l excitation and deexcitation (Draine 
Be |:toldi 1996| ). The H nucleon density was assumed 
to be nn = 10^ cm~^ (as expected for pressure bal- 
ance with the HII region), and the radiation field at 
1000 A w as assumed to b e x = 400 times stronger 



than the Habing (1968)| value (as expected for the 
~ 1.7 pc distance trom the exciting B0.5V star). 

The resulting model is shown in Figure ||, where 
it is assumed that cosO = 0.1, where is the an- 
gle between the line-of-sight and the normal to the 
model PDR. It is seen that the level populations 
are in fairly good agreement with observations, al- 
though the model overestimates the v = 1 and 2 
level populations, while underestimating the v = 3 
levels. For this model, T = 500 K at the point where 
y = 2n(H2)/nH = 0.2, and T = 150 K where y = 0.5. 
I I 

Timmermann et al. (1996) also compared predicted 
and observed fine structure line intensities with those 
predicted by the model. Agreement is generally good, 
allowing for uncertainties in abundances of Si and 
Fe: [CI]609/im is close to the observed value, while 



5. THERMAL BALANCE IN PDRS 



The obvious question now is: what temperatures do 
we expect the gas to have in PDRs, if we assume 
that the gas temperature is determined by a balance 
between heating and cooling? The principal heating 
processes in PDRs are: 

• Photoelectrons emitted from dust grains; 

• UV pumping, followed by collisional deexcita- 
tion of rovibrationally excited H2 ; 

• Photodissociation (kinetic energy of H+H); 

• H2 formation (kinetic energy, collisional deexci- 
tation); 

• Photoionization of H, C, etc. 



while the dominant cooling processes are: 

• fine-structure line emission from [CII]158/im, 
[OI]63,146/im, [SiII]35/im, [FeII]25/im,...; 

• H2 quadrupole emission; 

• CO rotational emission. 



The grain photoelectric heating rate is quite uncer- 
tain. Here we sho w results assuming the p hotoelec- 
tric heating rate of Bakes fc Tielens (1994)| . 



The cross sections for collisional excitation and deex- 
citation of H2 are also uncertain; for T^^o > 1000K 
we use the _H-H7 inelastic rates of [Martin fc Mandy 



1995 



and 



Mandy fc Martin (1995| 



while for 7 gas < 



000 K we extra polate the Ma,ndy fe Mart in rates us 
ing eq. (16) of [Praine fc Bertoldi (1996) , but using 
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Figure 5. Thermochemical profile of SI 40 PDR 
model, with temperatures determined by thermal bal- 
ance calculations. Lowel panel: gas temperature 
T, H2 fraction 2n(H2)/nH; and fractional ionization 
n(e)/nH. Upper panel: heating (solid lines) and cool- 
ing (broken lines) due to grain photoelectric heating 
(^^pe^^), the H2 molecule (including formation, colli- 
sional excitation and deexcitation, and photodissoci- 
ation), fine structure line emission ('fs."), and op- 
tical line emission (^'opt^'). Note that H2 has a net 
heating effect in much of the PDR. 



= 1000 K rather than 6 00 K as the transition tern - 
perature, as discussed by Draine & Bertoldi (1999). 



The fine-structure line emission is computed using 
standard collision rates, abundances characteristic of 
the gas toward COph (e.g., Si and Fe depleted by fac- 
tors of 40 and 250, respectively). Some of the lines - 
in particular, [01] 63 /im - can become optically thick. 
An "escape probability" ap proximation was used t o 
estimate local cooling rates ( [Draine fc Bertoldi 1999| ). 

Figure H shows the calculated thermochemical profile, 
as well as showing the contributions of different pro- 
cesses to heating and cooling. Close to the ionization 
front the fractional ionization is a few percent, and 
^gas ^ 4000 — 6000 K; in this region, optical line emis- 
sion (e.g., [SII]6716,6739) dominates the cooling, but 
is unimportant when Tgas ^ 4000 K. Emission from 
excited fine-structure levels dominates the cooling for 
50 < Tgas ^4000 K. 

H2 contributes net heating throughout much of the 
PDR. This heating arises mainly from collisional de- 
excitation of rovibrationally excited H2 following UV- 
pumping or H2 formation on grains, with a secondary 
contribution from the kinetic energy of newly- formed 
H2 and H+H following photodissociation. The other 
major heating process is photoelectric emission from 
dust grains. Both H2 heating and photoelectric heat- 
ing are important for A^h ^ lO^^cm"^; beyond this 
point H2 self-shielding has reduced the photodisso- 



Figure 6. Column densities N{v^ J) along the line- of - 
sight, divided by level degeneracies g{J). Solid sym- 
bols: observed values based on line intensities cor- 
rected for foreground extinction with Ak = 0.2mag^ 
and internal extinction within the PDR, inclined with 
with cosO = 0.1. Open symbols: model for SI40 PDR 
with X = 400 and nn = lO^cm"^^ with temperatures 
determined by thermal balance. 



elation rate (and therefore the H2 formation rate as 
well) to very low levels, and photoelectric heating 
dominates. 

Figure || shows the observed column densities of ex- 
cited states of H2 along the line of sight toward 
the S140 PDR, corrected for an assumed extinction 
Ak = 0.3 mag. A number of H2 line detections and 
upper limits were obtained subsequ ent to preparation 



of the Timmermann et al. (1996) paper on the S140 
PDR; Figure ^ shows column de nsities obtained from 



the most recent data reduction ( Bertoldi et al. 1998| ). 
Also shown are the column densities predicted by our 
model, which assumes an inclination of the PDR with 
0.1 - a factor of 10 "limb-brightening". 



cos 



There are significant differences between the model 
H2 level populations and those observed toward S140. 
The assumed factor of 10 limb-brightening leads to a 
good match for the vibrationally excited levels. How- 
ever, the V = column densities predicted by the 
model fall off too steeply - while the J = 2, 3 column 
densities in the model exceed the observed values by a 
factor ^ 2, the J = 3 — 7 column densities are too low 
by factors ~ 2 — 3. This presumably indicates that 
either the H2 inelastic cross sections or the thermal 
profile are incorrect. Note, however, that the model 
appears to be able to reproduce the observed popu- 
lation of H2 in J = 15, and thereby a high excitation 
temperature for the high— J levels. 

Somewhat better agreement can be obtained if the 
gas density and illuminating radiation field are both 
increased; in Figures ^ and g we show the result of in- 
creasing the illuminaung radiation field to x = 600, 



and raising the gas density to nn = 4 x lO^cm"^. 
With the increased x ^^^ nn, a good fit to the 
observed H2 line intensities is obtained for a hmb- 
brightening factor l/cos6> = 5. However, the as- 
sumed value of X is somewhat larger than expected 
from our knowledge of the exciting star, and nu 
is larger than expected given o ur estimates of the 
pressure in the ionization front ( [Timmermann et al. 
I9Q. It should also be noted that even with these 
values of x and nn, there are still noticeable discrep- 
ancies between model and observation. It is difficult 
to match the low observed line ratio between O-OS(l) 
and 0-0S(2). Both lines are in the same size ISO 
aperture, so beam dilution differences could not be 
the cause of the mismatch. The model also overpro- 
duces J = 15, although this line was observed after 
ISO^s helium boiloff, which made its flux calibration 
more difficult. Most other levels' column densities 
are well reproduced. 



erg s"VH 




0.01 



1018 1019 1020 

nodel 4e4_6e2_le4_3_6_5_10 iNnV^-iii J 



Figure 7. Same as Fig. 0^ but for nn = 4 x 10^ cm 
and X = 600. 



The high gas temperatures found in S140 are not 
unique. Early evidence for high gas temperatures 
in PDRs was obtain ed from II9 J = 3-1 and 4-2 
line observations by Parmar et al. (1991), who in- 



ferred Tgas = 1000 ± 2b[) K at a position .0244 pc 
from the ionization front in the Orion Bar, and 
^gas = 600 ±50 K another .0044 pc farther away from 
the ionization front. 



High gas temperatures near the H/H2 transition 
were inferred, based on ground- based spectra, for 



the reflection nebula NGC 2023 ( Draine fc Bertoldi 
19f6[). /go spectra of the reflection nebula NGC 70^ 
dBeirtoldi ^ ^^^^' ' ' ^ ' ^ ^ "^^ 
eTal 



i^oldi et al. 1999| ) and the Orion Bar ([Rosenthal 
/. 199^1 ), also show that some of the H2 must be 



in a region of gas temperatures T^ 



gas 



600 K. 
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Figure 8. Same as Fig. ^, hut for nn = 4 x 10^ cm~^ 
and X = 600; and with the PDR assumed to he in- 
clined with cosO = 0.2. 



WHAT NEEDS TO BE CHANGED IN THE 
MODELS? 



Where do we stand on modeling the H2 emission from 
PDRs? Considering that the emission arises from a 
region with a range of densities, one might take the 
point of view that we are doing pretty well - our 
constant density planar PDR models can approxi- 
mately reproduce the observed level populations us- 
ing ultraviolet pumping of H2 plus collisional exci- 
tation/deexcitation, with gas temperatures obtained 
from a balance of heating and cooling. Nevertheless, 
it seems to us that the discrepancies between our best 
models and the observations are unacceptably large. 
The discrepancies may be due in part to our under- 
lying assumptions of planar geometry and uniform 
density, but may also be due to poor approximations 
in our description of the local microphysics of heat- 
ing, cooling, chemistry, and level excitation. Below 
we indicate some of the areas where existing models 
may require modification. 



6.1. H2 inelastic collision cross sections 

The H-H2 collisional excitation/deexcitation rates 
used here for T^as ^ 600 K are those of [Martin fc 



Mandy (1995)| and [Mandy fc Martin (1995)| , which 
were obtained from quasi-classical trajectory calcu- 
lations. These rates are thought to be accurate for 
^gas ^ 600 K, but are likely to be increasingly inac- 
curate at low temp eratures, where quantum calcula - 
tions are required (t^lower 1997 , Forrey et al. 1997 ). 
Forrey et aL (1997)1 have recently reported accurate 
rates tor A J = ±2 transitions among the v = 0^ 

our extrap- 
results was 



J = - 5 levels, for 100 < T^^s < 1000 K: 



olation from the Mandy fc Martin (1995 



adjusted to give good agreement with these newer 
rates. 

While H-H2 inelastic rates are generally larger than 
for H2-H2 collisions, the latter are important when 
the H fraction becomes small. Qiiantum calcu lations 
have recently been reported by [Flower (1997)| . 



6.2. Grain Photoelectric Heating Rate 



The thermal models reported here em ploy a phc 
electric hea ting rate recommended by pakes fc Tie-| 



lens (1994) for heating by 6-13.6 eV photons: 



^YL^uv^uv^ ^ 



6 X 10~^^cm^ 



auv — dust EUV absorption/H 
Uuv = EUV energy density = X x 0.033 eVcm~^ 
e = photoelectric heating efficiency 

= .049(1. + x/1925)"^-^^ + 
.037T4°-'^(l. + x/5000)-^ 
X ^ x{T/Kf-'/{nJcra-^) 



Bakes & Tielens assumed a "MRN" power-law size 
distribution down to ultrasmall grains containing 
only ^ 30 carbon atoms (i.e., PAH molecules). Since 
we have only an imprecise understanding of the com- 
position and size distribution of the grains, this heat- 
ing rate is necessarily uncertain. Under the condi- 
tions of interest (e.g., the point where 2n(H2)/nH = 
0.1 in Figure pi), we have x ^ 400, ng ~ 1.5 cm~^, 
and T ^ 350 k, so that x ^ 5000, and the heat- 
ing efficiency e ^ 0.021, which increases further into 
the cloud to e ?^ 0.049. Since photoelectric emission 
from grains is a dominant heating process in PDRs, 
its further study would be of value. 



6.3. Dust-to-Gas Ratio 



The dust grains play a crucial role in photodissocia- 
tion regions. In particular, they determine the rate 
for formation of H2 from H, and, as seen above, pho- 
toelectric emission from dust is a major heat source. 
Common practice is to adopt a "standard" dust-to- 
gas ratio and wavele ngth-dependent dust opacitv in 



the PDR. However, |Weingartner fc Draine (1999) 
have pointed out that the anisotropic radiation in a 
PDR can drive the grains through the gas at drift ve- 
locities which may be comparable to the flow velocity 
of the gas relative to the dissociation front, thereby 
leading to an increase in the dust-to-gas ratio in the 
PDR, at least for those grain compositions and sizes 
for which the grain-gas drift velocities are large. In- 
creased dust concentrations in part of the PDR could 
alter the thermochemical properties of the PDR by 
increasing the H2 formation rate (and therefore the 
local rate of UV pumping and photodissociation of 
H2) and by increasing the grain photoelectric heat- 
ing rate. As a result, we can expect that the gas 
temperature will be increased near the photodissoci- 
ation front. 



The anisotropic radiation incident on a dust grain re- 
sults in a force which consists of the "ordinary" radia- 
tion pressure force due to absorption and scattering, 
plus additional forces due to anisotropic photoelec- 
tron emission (the "illuminated" side of the grain will 
emit more photoelectrons per second than the "dark" 
side) and anisotropic photodesorption (adsorbed H 
or H2 can be photodesorbed, with a larger rate on 
the "illuminated" side of the grain). The recoil force 
from photoelectric emission and photodesorption can 
be several times larger than the "ordin ary" radiation 
pressure (Weingartner & Draine 1999). 



Because the grains are charged, drift across the mag- 
netic field lines will be inhibited. It will be of great 
interest to see the extent to which these effects will 
change the dust-to-gas ratio, and the grain size dis- 
tribution, in the PDR, and how this affects the ther- 
mochemical profile of the PDR. 



7. SUMMARY 



ISO has opened a new window on photodissociation 
fronts, by allowing us to measure the populations of 
the rotationally-excited levels of the v = state of 
H2, thereby placing strong constraints on the gas 
temperature. The S140 PDR provides an excellent 
illustration of this capability. 

The observed line emission from rotationally- 
excited levels of H2 indicates that gas temperatures 
^gas ^ 500 K prevail in part of the PDR where ap- 
preciable H2 is present. Attempts to construct theo- 
retical models of the PDR tend to fall short in terms 
of the populations of levels J = 3 — 7, 

We discuss the uncertainties in the dominant heating 
and cooling mechanisms; the H2 inelastic collision 
rates still have significant uncertainties, and the grain 
photoelectric heating rate is yet not well-established. 
In addition, it appears possible that the dust-to-gas 
ratio in PDRs such as S140 may deviate from the 
"standard" value since the anisotropic radiation field 
can drive the grains through the gas with significant 
drift speeds. 
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